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ABSTRACT 

We present results of a variability search in the intermediate-age open cluster 
NGC 6866 from 29 nights over two observing seasons. We find 28 periodic variables, of 
which 19 are newly identified. The periods of these variables, which have V magnitudes 
from 11.5 to 19.3 mag, range from ~ 48 min to 37 d. We detected several J-Scuti stars, 
some of which are of high amplitude, as well as 7-Doradus, rotational variables and 
eclipsing binaries. In order to study the physical properties of the cluster, we obtained 
U BVRI photometry of all the stars on a good photometric night. The radial distribu- 
tion of the stellar surface density shows that the cluster has a radial extent of about 7 
arcmin (~ 3 pc) with a peak density of 5.7 ± 0.7 stars/arcmin 2 at the cluster center. 
The colour-colour diagram indicates a reddening of E{B — V) = 0.10 mag towards 
NGC 6866. A distance of ~1.47 kpc and an age of ~ 630 Myr is estimated from the 
colour-magnitude diagram using the theoretical isochrones of solar metallicity. 

Key words: open cluster:individual:NGC 6866-stars: variables : general- 
techniques :photometric 



1 INTRODUCTION 

The study of Galactic open clusters is important for under- 
standing the history of star formation and the nature of the 
parent star clusters. These systems are used to test stellar 
models and are vital for our understanding of stellar evolu- 
tion. By comparing the colour-magnitude diagram of a star 
cluster with theoretical evolutionary models, its age may 
be estimated and information on the evolution of stars of 
nearly the same age and chemical composition could be ob- 
tained. However, uncertainties in various parameters such as 
reddening, distance and chemical composition compromise 
our efforts to measure the age and other parameters. Some 
cluster members also show various kinds of variability at dif- 
ferent stages of their evolution. For example, a wide variety 
of pulsating variable stars such as 8 Scuti and /3 Cephei stars 
are found in open star clusters. Observations of these vari- 
ables can be used to infer many important stellar parameters 
such as their masses, radii and luminosities. These param- 
eters are required for an investigation of the relationships 
between rotation, stellar activity, age and masses (Messina 
et al. 2008, 2010) and impose constraints on stellar pulsation 
models. Recent work on cluster variable stars can be found 
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in Hargis & Sandquist (2005), Parihar et al. (2009), Marchi 
et al. (2010), Saesen et al. (2010), among others. 

At ARIES, Nainital, we have started a long-term ob- 
servational survey of variable stars in unstudied or poorly 
studied young and intermediate age open clusters. The pur- 
pose of this survey is to search and characterize the vari- 
able star content of these clusters and to determine their 
fundamental parameters. Young clusters allow the study of 
pre-main sequence (PMS) stars and the effects of variabil- 
ity on the spread in cluster ages. Intermediate age clusters 
can be used to study short-period variables such as S Scuti 
and 7 Dor stars. Intensive un-interrupted observations over 
a few hours, as well as extensive observations spanning sev- 
eral months and years, are planned for many such clusters. 

In this paper, we present a study of the intermediate age 
open cluster NGC 6866 (RA = 20:03:55, DEC = +44:09:30; 
I = 79°. 560, b = +6°. 839). The cluster is in the field ob- 
served by the Kepler satellite, whose aim is to search for 
earth-like planets by monitoring the light variations of more 
than 150 000 stars in a 105 square degree field situated in 
the Cygnus-Lyra region. The data obtained in this mission 
are also used for the asteroseismic study of pulsating vari- 
ables and to investigate stellar characteristics with micro- 
magnitude precision and uninterrupted time coverage (see, 
Basri et al. 2010). The cluster NGC 6866 has not been ex- 
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tensively surveyed for variability. The Hipparcos and Tycho 
catalogues contain a few suspected variables in the cluster. 
Molenda-Zakowicz et al. (2009) obtained time-series CCD 
photometric observations of this cluster for 14 nights over a 
period of about 3 months. They found 19 variables in the 
field of the cluster, of which 12 seem to be periodic variables. 
However, due to their limited time coverage, the periods and 
characteristics of these variables are somewhat uncertain. 

Details of the observations and data reduction are pre- 
sented in §2. The photometric calibration and cluster pa- 
rameters are studied in §3. In §4 we discuss identification of 
variable stars and their characteristics. This is followed by 
a discussion and the conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

Time-series observations of stars in NGC 6866 were obtained 
in the Johnson V and Cousin I bands using the 1.04-m Sam- 
purnanand telescope at Manora Peak, Nainital. We used a 
2k x 2k CCD camera covering a field of view of about 13 
x 13 arcmin. The readout noise and gain of the CCD are 
5.3 e- and 10 e-/ADU respectively. The cluster was observed 
for 29 nights between 2008 September 26 and 2011 Jan- 
uary 10 over two observing seasons. Priority was given to 
F-band observations in order to increase the time sampling 
in this waveband. A total of 768 frames in the V band and 
50 frames in the I band were accumulated. Bias and twi- 
light flat frames were taken on a regular basis. To improve 
the signal-to-noise ratio, all observations were taken in 2 x 2 
binning mode (~ 0.74 x 0.74 arcmin). The exposure time 
for each frame ranged from 30 sec to 200 sec depending on 
sky conditions and time constraints. Typical seeing was ~ 2 
arcsec. In order to detect short-period variables we observed 
NGC 6866 continuously for more than 4 hours on 3 different 
nights with an observing cadence of about 2.2 min. A brief 
log of the observations is given in Table 1. 

The basic steps of image processing, which include bias 
subtraction, flat fielding, and cosmic ray removal, were per- 
formed using IRAF. During the observations, it was not pos- 
sible to keep all the stars at the same pixel position on ev- 
ery night. The IRAF tasks GEOMAP and GEOTRAN were 
used to align all the images with respect to a reference frame 
which was chosen so that the cluster center falls approxi- 
mately at the center of the observed field. In most cases, we 
achieved a pixel accuracy of ?a 0.1 arcsec. 



3 PHOTOMETRIC STUDY OF THE CLUSTER 
NGC 6866 

3.1 Photometric Calibration 

Photometry was performed using DAOPHOT II profile fit- 
ting software (Stetson 1987). In order to perform consistent 
photometry from night to night on the aligned images, we 
made a master list of 2809 stars from 13 frames which have 
the best seeing and coverage of the target field. It should be 
noted that not all stars could be measured on every frame 
owing to different exposure times and sky conditions. 

To determine the extinction coefficients, colour terms 
and zero points, we used Landolt's (1992) standard fields 
SA95 and PG 0231+051. These were observed on the night 



Table 1. Log of the CCD observations. 
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of 2010 November 30 in photometric sky conditions. Expo- 
sure times of 300, 300, 200, 100, and 60-sec were used for U, 
B, V, R and I respectively. We acquired two frames each 
of NGC 6866 and PG 0231+051 and four frames of SA95. 
A total of fifteen standard stars with 12.68 < V < 16.28 
and —0.53 < (V — I) < 1.95 were observed during this 
night with airmasses ranging from 1.11 to 2.04. Nine of 
these stars are in SA95 and six in PG 0231+051. We used 
profile-fitting photometry to extract the magnitudes using 
the program DA OGR OW for obtaining the curve-of-growth 
corrections. A linear least square regression was used to 
the standard system to derive the following transformation 
equations: 

u = U + 8.13 ± 0.01 - (0.06 ± 0.01) (17 - B) + (0.59 ± 0.04)A 
b = B + 5.84 ± 0.02 - (0.02 ± 0.02) (S -V) + (0.27 ± 0.03)X 
v = V + 5.40 ± 0.01 - (0.02 ± 0.01)(B -V) + (0.17 ± 0.02)X 
r = R + 5.22 ± 0.01 - (0.04 ± 0.02) (V - R) + (0.09 ± 0.02)X 
i = I + 5.63 ± 0.02 + (0.01 ± 0.01) (Ji - I) + (0.05 ± 0.02)A 

where u,b,v,r and i are the aperture instrumental mag- 
nitudes and U, B, V, R and I are the standard magnitudes 
and X is the airmass. Zero-point and colour-coefficient er- 
rors are ~ 0.01 mag. These equations were used to gen- 
erate secondary standard stars in the target field observed 
on the same night. To standardize the data on remaining 
nights, differential photometry was performed using these 
secondary stars. We used a linear fit between the stan- 
dard and instrumental magnitudes on each night, assuming 



Photometric Study and Detection of Variable 



Stars in the Open Clusters - I : NGC 6866 3 




12 14 16 R 18 20 22 



b 




12 14 16 \; 18 20 22 



b 




12 14 16 g 18 20 22 




Figure 1. The photometric error of stars as a function of bright- 
ness in UBVRI bands. 

that most of the stars are non-variables. We rejected stars 
which deviated by more than 3 — a deviations. In Fig. 1, we 
show the standard deviations as a function of magnitudes in 
UBVRI bands. The photometric errors were computed on 
the night of standardization. The typical standard deviation 
is less than 0.05 mag for stars brighter than ~ 19 mag in V 
band. 

3.2 Comparison with previous photometry 

NGC 6866 has not been intensively observed in the past. 
UBV photoelectric and photographic photometry were ob- 
tained by Hoag et al. (1961), Johnson et al. (1961), and 
Barkhatova and Zakharova (1970), mostly for bright stars. 
Recently, Molenda-Zakowicz et al. (2009; hereafter MOL09) 
obtained BVIc CCD photometry in a 12.8 x 11.7 arcmin 
field centred on NGC 6866 using the 60-cm telescope at 
the Bialkow Observatory of the Wroclaw University, 
Poland. They observed 552 stars and derived calibrated B 
and V photometry. We did not compare our photometry 
with the older photoelectric and photographic work as these 
contain relatively large errors and are mostly confined to the 
brighter stars. We do, however, compare our results with 
MOL09 which, in turn, have been compared with older pho- 
toelectric and photographic measurements. We have cross- 
identified our stars with those of MOL09 and find 511 stars 
in common. Fig. 2 shows the difference, A, in V and (B — V) 
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Figure 2. A comparison of our V and (B — V) observations with 
those of MOL09. The solid line represents a zero magnitude dif- 
ference while the dashed line in the lower panel is a least-squares 
linear fit. 



between our photometry and MOL09. A linear fit between 
our V magnitudes and those of MOL09 shows a magnitude 
dependence which can be given by the following equation 

Vpresent ~ V Mo l = (0.008 ± 0.001) X Vpresent ~ (0.13 ± 0.02) 

Our V magnitudes are noticeably brighter than 
MOL09 for the brightest stars, particularly towards 
V < 13 mag. For faint stars, the V magnitudes and (B — V) 
colours between two photometry seem to be in agreement 
within their internal photometric errors. 



3.3 Astrometry 

To convert the pixel coordinates (X,Y) into celestial coordi- 
nates (cV20oo, 52ooo ), a linear astrometric solution was derived 
by matching about 134 bright stars in common between our 
catalogue and the 2MASS catalogue. We used the CCMAP 
and CCTRAN tasks inside IRAF for this purpose. Using 
these equations, we converted the (X,Y) position into J2000 
celestial coordinates for all the stars in our catalogue. The 
typical rms scatter, as derived from the difference between 
calibrated and catalogued coordinates, is about 0.1 arcsec. 

In Table 2, we present UBVRI photometric measure- 
ments for 2473 stars identified in the field of the cluster 
NGC 6866. This is only a sample listing; the full table is 
available in electronic form. The table contains the star ID, 
its celestial coordinates for J2000 and UBVRI magnitudes 
and their corresponding errors. It should be noted that we 
were able to determine the U magnitudes for only about 
half of these stars owing to the low quantum efficiency of 
the CCD in the U band. 
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Table 2. UBVRI photometric catalogue of 2473 stars detected in the field of cluster NGC 6866. The error in magnitudes indicates 
the internal photometric error in the measurement. Table is sorted in the order of increasing V magnitude. Last 3 columns give the 
membership probabilities based on the spatial distribution, their position in the CMD and proper motion. 



ID 



RA 



DEC 



U 



B 



V 



R 



1 



sp 



' P h 



! pm 



l 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 



20:04: 
20:03: 
20:04: 
20:03: 
20:04: 
20:03: 
20:03: 
20:03: 
20:03: 
20:04: 
20:03: 
20:03: 
20:03: 
20:03: 
20:03: 
20:03: 
20:03: 
20:03: 
20:04: 
20:04: 
20:03: 
20:03: 
20:03: 
20:03: 
20:04: 
20:04: 
20:03: 
20:04: 
20:03: 
20:04: 



06.21 
52.44 
09.32 
57.62 
02.85 
52.57 
57.35 
40.66 
52.57 
22.90 
54.70 
22.82 
55.17 
51.69 
57.70 
26.12 
55.18 
53.15 
27.27 
25.52 
23.96 
56.83 
39.17 
49.52 
06.14 
18.38 
47.13 
14.16 
28.54 
03.97 



44:04:55.0 
44:11:31.0 
44:04:16.2 
44:08:37.5 
44:11:55.4 
44:11:33.1 
44:09:33.5 
44:11:43.5 
44:15:46.9 
44:10:31.2 
44:10:52.9 
44:15:50.3 
44:08:24.2 
44:10:19.1 
44:06:27.6 
44:10:05.3 
44:10:44.0 
44:09:20.4 
44:10:48.8 
44:10:16.2 
44:10:46.8 
44:10:32.0 
44:14:19.1 
44:10:50.7 
44:12:54.5 
44:09:52.2 
44:09:25.7 
44:12:32.9 
44:09:42.1 
44:10:20.5 



11.70±0.01 
10.92±0.01 
10.98±0.01 
11.27±0.01 
11.52±0.01 
11.17±0.01 
11.59±0.01 
10.92±0.01 
14.90±0.01 
11.77±0.01 
11.80±0.01 
13.29±0.01 
13.26±0.01 
12.06±0.01 
13.38±0.01 
12.11±0.01 
12.29±0.01 
12.33±0.01 
14.03±0.01 
12.45±0.01 
14.49±0.01 
12.41±0.01 
12.50±0.01 
12.39±0.01 
14.22±0.01 
12.52±0.01 
12.65±0.01 
12.87±0.01 
14.36±0.01 
12.74±0.01 



ll.01i0.01 
10.75±0.01 
10.95±0.02 
11.08±0.01 
11.31±0.01 
11.15±0.02 
11.43±0.01 
11.23±0.01 
12.90±0.01 
11.56±0.01 
11.64±0.01 
12.52±0.01 
12.64±0.01 
11.94±0.01 
12.73±0.01 
12.06±0.01 
12.14±0.01 
12.36±0.01 
13.10±0.01 
12.30±0.01 
13.27±0.01 
12.27±0.01 
12.30±0.01 
12.26±0.01 
13.24±0.01 
12.36±0.01 
12.52±0.01 
12.81±0.01 
13.39±0.01 
12.63±0.01 



10.03±0.07 
10.57±0.04 
10.62±0.01 
10.75±0.03 
11.08±0.01 
11.09±0.01 
11.19±0.01 
11.22±0.01 
11.29±0.03 
11.35±0.01 
11.37±0.01 
11.53±0.03 
11.60±0.01 
11.63±0.01 
11.69±0.01 
11.72±0.01 
11.86±0.01 
11.95±0.01 
11.97±0.01 
11.98±0.01 
11.99±0.01 
12.02±0.01 
12.02±0.01 
12.04±0.01 
12.07±0.01 
12.08±0.01 
12.20±0.01 
12.21±0.01 
12.28±0.01 
12.31±0.01 



9.54±0.08 
10.40±0.03 
10.44±0.01 
10.59±0.02 
10.94±0.01 
11.01±0.03 
11.08±0.01 
11.20±0.01 
10.44±0.02 
11.23±0.01 
11.21±0.01 
10.98±0.03 
11.09±0.01 
11.48±0.01 
11.19±0.01 
11.51±0.01 
11.74±0.01 
11.74±0.01 
11.39±0.01 
11.81±0.01 
11.33±0.01 
11.91±0.01 
11.85±0.01 
11.94±0.01 
11.46±0.01 
11.91±0.01 
12.02±0.01 
11.87±0.01 
11.72±0.01 
12.14±0.01 



9.13±0.03 
10.26±0.03 
10.26±0.01 
10.40±0.01 
10.80±0.01 
ll.01i0.01 
10.97±0.01 
11.22±0.01 
9.55±0.03 
11.11±0.02 
11.06±0.01 
10.53±0.06 
10.58±0.01 
11.31±0.01 
10.69±0.01 
11.32±0.04 
11.61±0.01 
11.51±0.02 
10.81±0.01 
11.60±0.01 
10.69±0.03 
11.78±0.01 
11.70±0.02 
11.82±0.02 
10.88±0.01 
11.73±0.01 
11.82±0.01 
11.54±0.01 
11.19±0.02 
11.95±0.01 



0.26 
0.72 
0.14 
0.83 
0.60 
0.72 
0.91 
0.54 
0.12 
0.25 
0.82 
0.00 
0.83 
0.88 
0.54 
0.28 
0.84 
0.96 
0.13 
0.19 
0.21 
0.85 
0.23 
0.79 
0.43 
0.37 
0.82 
0.33 
0.35 
0.72 



1.00 
1.00 
1.00 
1.00 
1.00 
1.00 
1.00 
1.00 
0.00 
1.00 
1.00 
1.00 
1.00 
1.00 
1.00 
1.00 
1.00 
1.00 
0.00 
1.00 
0.00 
1.00 
1.00 
1.00 
0.00 
1.00 
1.00 
1.00 
0.00 
1.00 



0.00 
0.76 
0.00 
0.91 
0.82 
0.76 
0.89 
0.94 
0.65 
0.64 
0.54 
0.88 
0.73 
0.68 
0.26 
0.66 
0.53 
0.13 
0.65 
0.16 
0.72 
0.59 
0.93 
0.74 
0.02 
0.40 
0.80 
0.00 
0.02 
0.83 



2472 
2473 



20 03 33.64 
20:03:51.71 



44 13 43.2 
44:10:37.0 



21.85±0.31 
22.02±0.49 



21.24±0.30 
21.50±0.38 



20.19±0.20 
20.51±0.26 



18.99±0.32 
19.53±0.19 



0.21 
0.84 



0.08 
0.22 



3.4 Cluster parameters 

3.4-1 Radius 

The spatial structure of an open cluster, as derived from the 
stellar density distribution obtained by star counts, is diffi- 
cult to estimate because of the irregular shape of the cluster. 
However, we can obtain some information about the core 
and the approximate cluster radius. To determine the radial 
density distribution (RDP), we divided the cluster into a 
number of concentric annular regions centred on the cluster. 
The center of the cluster was defined by the pixel coordinate 
(483, 519) in our reference CCD frame. The cluster radius 
is defined as the radius where the surface density becomes 
approximately equal to the average density of the surround- 
ing field. To reduce field star contamination, we only used 
stars brighter than 19 mag in V in determining the stellar 
surface density. The spacing of the annular regions was set 
to 40 pixels (~ 30 arcsec). This was chosen in such a way 
that each zone contains a statistically significant number of 
stars. 

We counted the number of stars in each annular region. 
The number density in i th zone is then Pi = ^r, where iV» 
and Ai are the number of stars and the area of the i th zone 
respectively. The error bars are derived assuming that the 



number of stars in a zone follows Poisson statistics. The 
radial variation of the derived density distribution is shown 
in Fig. 3. 

Kaluzny (1992) models the projected stellar density dis- 
tribution as 

P(r) = Pf + T+T%F 

where po is the density at the cluster center, pf is the field 
density and r c is the core radius of cluster. The core radius 
is defined as the radius where the density, p(r), is half of 
the central density, po. A fit to this function is shown by the 
solid line in Fig 3. Using \ 2 minimization, we determined 
the core radius and central density to be 2.0 ± 0.5 arcmin 
and 5.7 ± 0.7 stars/arcmin 2 respectively. Even though our 
CCD frames did not encompass a large area, we are able to 
determine that the surface density is approximately equal 
to the density of the surrounding field stars at a radius of 
about 7 arcmin, which we take as the cluster radius. 

3.4.2 Reddening 

The reddening, E(B — V), in the cluster region can be de- 
termined using the (U — B) vs (B — V) two-colour diagram 
(TCD). Because only half the stars have U measurements, 
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Figure 3. The stellar radial density distribution for stars brighter 
than 19.0 mag in the field of NGC 6866. The solid line represents 
the King profile and the horizontal dashed line shows the average 
field star density. 

and also to avoid contamination by field stars, we used only 
stars with V < 17 mag and au, ob, <?v smaller than 0.05 
mag. The resulting TCD is shown in Fig. 4. The observed 
cluster sequence is clearly seen, as is the turn off which oc- 
curs at about spectral type AO. The most probable cluster 
members with V < 17.0 mag and (B — V) < 0.9 mag are 
shown by a different symbol. 

We adopted the slope of the reddening vector ^(b-v) = 
0.72. We then fitted the intrinsic zero-age main sequence 
(ZAMS) with solar metallicity, as given by Schmidt-Kaler 
(1982), to the observed main-sequence (MS) stars. This was 
done by shifting E(B — V) and E(U — B) along the reddening 
vector (shown in Fig. 4 by the dashed line). The best match 
results in a mean reddening of E(B — V) — 0.10 mag for 
NGC 6866. 

3.4-3 Distance and Age 

The identification of the cluster main sequence on the colour- 
magnitude diagrams (CMDs) allows a model-dependent 
mass, radius, and distance for each star to be determined. 
These are determined by fitting isochrones assuming that 
the star is on the main sequence. To derive these parame- 
ters, we constructed calibrated (B — V), (V—R), and (V — I) 
vs V diagrams of NGC 6866 using our data (Fig. 5). A well- 
defined cluster main sequence is clearly seen in all the CMDs. 
However, they are contaminated by the typical red field star 
population. In order to determine the age and distance of 
the cluster, theoretical isochrones of Girardi et al. (2002) for 
solar metallicity are over-plotted on the CMDs (solid lines in 
Fig. 5). Girardi isochrones are visually fitted by varying the 
distance modulus and age simultaneously while keeping the 
reddening fixed at E(B — V) = 0.10 mag (as estimated ear- 
lier using TCD). We have assumed that the total extinction 
Av = 3.1 x E(B — V). The best-fit model of the cluster gives 
a reddening-free distance modulus (Vo — Mv) = 10.84 mag 
and an age of log(t) = 8.8 yrs. 

The physical parameters determined for the cluster 
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Figure 4. The (U — B) versus (B — V) diagram for stars in 
NGC 6866. The larger filled circles are stars which lie on or close 
to the main sequence in the colour-magnitude diagram and which 
have proper motions suggesting more than a 60 percent probabil- 
ity of being cluster members. The dashed straight line represents 
the slope (0.72) and direction of the reddening vector. The solid 
curve represents the zero-age main sequence from Schmidt-Kaler 
(1982). 



Table 3. A comparison of the fundamental parameters of 
NGC 6866 from our study and that of Kharchenko et al. (2005) 



Cluster parameter Present study Kharchenko et al. (2005) 



Rcore (arcmin) 


2.0 


4.2 


Rduster (arcmin) 


7.0 


8.4 


Mean E{B - V) 


0.10 


0.17 


V -M v 


11.15 


11.33 


log(Age/yr) 


8.8 


8.68 


Distance (kpc) 


1.47 


1.45 



NGC 6866 in the present study are summarized in Table 3. 
We also compare our results with those given by Kharchenko 
et al. (2005). In general, our values are in good agreement. 
However, we obtained a significantly smaller core radius and 
reddening than Kharchenko et al. (2005). 

4 VARIABLE STARS IN THE CLUSTER 
NGC 6866 

4.1 Selection criteria to identify variables 

After excluding observations taken under poor sky condi- 
tions, we considered only 718 photometric observations in 
V and 47 in I bands for the variability search. In order to 
search variable stars in the target field, the time-series V 
band magnitude of all the stars are passed through the pe- 
riodicity analysis explained in the following section. A star 
was analysed for variability only if it met the following cri- 
teria. 

(i) The star should not be within 10 pixels of the edge of 
any frame. 
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Figure 5. The colour-magnitude diagrams, (B — V), (V — R) and 
(V - I) as a function of V, for stars in the field of NGC 6866. The 
solid line is the Girardi et al. (2002) isochrone (solar metallicity) 
for a cluster age of log(i) = 8.8. A distance modulus of V — My 
= 11.15 and reddening E(B — V) = 0.10 mag was used. 



(ii) The magnitude of the star should have a standard 
deviation less than three times the mean error in the corre- 
sponding magnitude bin. 

(iii) The stellar magnitude needs to be measured in at 
least 100 frames. 



4.2 Periodicity analysis 

We used the Lomb-Scarglc periodogram (Lomb 1976, Scar- 
gle 1982) to estimate the period of a variable star. This 
method computes the Fourier power spectrum by fitting sine 
and cosine terms over a large number of frequencies in the 
given frequency range and is applicable to unevenly sampled 
data. The period was derived using only the V-band data 
because of the better sampling in this filter. In some cases, 
a few points were removed because either they were devi- 
ated by more than 0.2 mag in two sequential observations 
or strongly deviated from the mean magnitude. To search 
periodic variables, we have given a range of periods between 
0.01 days to 100 days as our total time length was about 100 
nights in 2010. The large range of period search was chosen 
because Lomb-Scargle method uses a period search in the 
logarithmic time scale hence shallower at the longer peri- 
ods. However, we have not considered any period beyond 
50 days, which is about half of the total observing span. 
We noticed that many spurious variables were found with 
periods which are harmonics or daily aliases of each other. 
These stars were rejected. The only stars lie in the magni- 
tude range V < 19.5 were considered for variability search 
as photometric magnitudes have large errors towards fainter 
end. Further, we have not considered 2008 observing data 
for the brighter stars having V < 12.5 mag which shows 
large uncertainty in the photometric measurements. 

After these selection criteria, we visually inspected the 
light curves of the remaining stars for periodic variation and 
found 28 variable stars. A phased light curve was derived for 
each star using an estimated period. A phase-folded light 
curve, binned in intervals of 0.01 phase, was constructed by 
taking the average of the magnitudes in these phase bins. 
It was noticed that some periodic variables, particularly in 
binary stars, the period-finding routine gives a 'best-fit' pe- 
riod that actually corresponds to half the orbital period. To 
detect such variables, we inspected the light phased at twice 
the period given by the period-finding algorithm. For exam- 
ple, variable ID 248 gives a period of 0.217-day using the 
Lomb-Scargle algorithm, but by visual inspection we found 
that the true period is 0.434-d because the light curve has 
two unequal minima. In this way, we found seven stars which 
were better represented with twice the period given by the 
Lomb-Scargle periodogram. 

Among the 28 periodic variables identified in the 
present study, 19 are newly detected. Table 4 lists the iden- 
tification number of variable stars, celestial coordinates, pe- 
riod, and its membership probabilities. We also give the 
intensity-averaged mean magnitude and the amplitude of 
brightness variation in V band. The amplitude was esti- 
mated as the difference between the median values of the 
upper and lower 15% magnitude values of the light curve 
(cf., Herbst et al. 2002; Messina et al. 2010). In the last two 
columns, we give a possible class of variability from our anal- 
ysis and corresponding identification in MOL09 and period 
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Figure 6. V band phase light curve for the 28 variable stars identified in the present study. Phase is plotted twice and shifted in such 
way that the light minimum falls near to zero phase. 
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Table 4. Details of the 28 periodic variables identified in the present study. The variable stars identification are sorted in the order 
of increasing IDs. The columns give RA, DEC, period, V band mean magnitude, (B — V) colour, amplitude of variation in V band, 
membership probabilities based on the spatial distribution, their position in the CMD and proper motion, cluster membership, and their 
possible variability class. Last column gives the characterization, identification number and period, if the star is already identified in 
MOL09. 



Star 


RA 


DEC 


Period 


< V > 




AV 


Membership Probability 


Classification based on 


ID 


(J2UUU J 


(J2UUU J 


(days) 


(mag) 


(mag) 


(mag) 


Psp 


Pph 


Ppp, 


Status 


present study 


MOL09 (ID, P) 


0016 


20:03:26.12 


44:10:05.3 


0.465333 


11.747 


0.34 


0.064 


0.28 


1.00 


0.66 


likely 


binary? 




0020 


20:04:25.52 


44:10:16.2 


0.380373 


12.003 


0.32 


0.032 


0.19 


1.00 


0.16 


unlikely 


binary? 




0027 


20:03:47.13 


44:09:25.7 


0.143616 


12.222 


0.32 


0.039 


0.82 


1.00 


0.80 


member 


S Set 


S Set (V3, 0.106±0.02) 


0036 


20:03:42.47 


44:10:06.4 


0.214684 


12.623 


0.32 


0.035 


0.70 


0.47 


0.84 


likely 


5 Set 




0039 


20:03:54.18 


44:06:46.0 


0.836120 


12.677 


0.25 


0.035 


0.60 


1.00 


0.81 


likely 


7 Dor 


7 Dor (V12, 0.707) 


0047 


20:04:11.20 


44:05:33.3 


0.275710 


12.995 


0.31 


0.041 


0.27 


1.00 


0.55 


unlikely 


? 


S Set (VI, 0.066) 


0058 


20:03:31.61 


44:07:59.8 


0.911577 


13.246 


0.97 


0.052 


0.38 


0.00 


0.00 


field 


7 Dor 




0074 


20:03:34.93 


44:14:50.1 


0.321750 


13.469 


0.43 


0.055 


0.11 


1.00 


0.72 


likely 


Ell 


W UMa (V5, 0.321) 


0081 


20:03:27.93 


44:09:19.1 


1.239157 


13.548 


0.43 


0.058 


0.33 


1.00 


0.92 


likely 


7 Dor 


Irr (V13) 


0094 


20:03:59.34 


44:10:25.8 


0.740741 


13.878 


0.43 


0.087 


0.82 


1.00 


0.77 


member 


7 Dor 


7 Dor (Vll, 0.805) 


0158 


20:03:40.87 


44:09:40.0 


2.285714 


14.902 


0.68 


0.052 


0.66 


1.00 


0.75 


member 


EA 




0191 


20:03:33.48 


44:13:53.4 


1.090513 


15.285 


0.67 


0.085 


0.19 


1.00 


0.67 


likely 


binary? 




0221 


20:03:53.35 


44:04:04.2 


0.476872 


15.564 


0.83 


0.079 


0.21 


1.00 


0.86 


likely 


PV 




0231 


20:03:48.28 


44:10:56.6 


0.327547 


15.615 


0.70 


0.052 


0.76 


1.00 


0.81 


member 


PV 




0239 


20:04:19.00 


44:07:05.8 


37.037037 


15.63 


2.09 


0.26 


0.27 


0.00 


0.81 


unlikely 


semi-reg 


Irr (V18) 


0248 


20:03:38.79 


44:04:53.0 


0.437446 


15.66 


0.82 


0.21 


0.22 


1.00 


0.74 


likely 


Ell 


Ell (V9, 0.434) 


0253 


20:04:20.84 


44:10:03.9 


7.407407 


15.73 


1.18 


0.06 


0.31 


0.00 


0.32 


field 


Rot 


Irr (V14) 


0332 


20:03:43.64 


44:05:19.7 


11.494253 


16.19 


1.44 


0.14 


0.33 


0.00 


0.18 


field 


semi-reg 




0349 


20:04:11.78 


44:13:06.9 


12.048193 


16.31 


1.41 


0.10 


0.32 


0.00 


0.97 


unlikely 


semi-reg 




0444 


20:04:21.25 


44:15:40.6 


16.260162 


16.82 


1.47 


0.31 


0.00 


0.00 


0.00 


field 


EB 




0487 


20:03:49.82 


44:11:08.5 


0.415110 


17.21 


0.93 


0.34 


0.76 


1.00 


0.85 


member 


W UMa 


W UMa (V7, 0.415) 


0494 


20:04:00.17 


44:14:03.2 


0.366704 


17.26 


1.03 


0.44 


0.34 


1.00 


0.00 


unlikely 


W UMa 


W UMa (V6, 0.366) 


1077 


20:04:13.87 


44:03:45.8 


0.033559 


18.58 


0.66 


0.34 


0.02 


0.00 


0.96 


unlikely 


HADS 




1088 


20:03:56.20 


44:12:49.9 


0.184775 


18.67 


1.02 


0.28 


0.53 


0.04 


0.03 


field 


HADS 




1274 


20:04:26.66 


44:05:35.9 


0.462428 


19.05 


1.42 


0.46 


0.00 


1.00 




unlikely 


W UMa 


W UMa (V4, 0.262) 


1292 


20:03:41.23 


44:12:17.9 


0.295247 


18.95 


1.72 


0.37 


0.49 


0.30 


0.98 


unlikely 


HADS 




1421 


20:03:41.08 


44:08:47.4 


0.041263 


19.15 


0.87 


0.30 


0.65 


0.00 


1.00 


unlikely 


HADS 




1583 


20:03:58.70 


44:11:33.5 


0.082055 


19.26 


0.99 


0.32 


0.70 


0.00 


0.73 


likely 


HADS 





if the star has already been identified. The V band phased 
light curves of these variables are shown in Fig. 6. 

Having the lesser number of data points, smaller ampli- 
tude of variation and larger photometric error in I band in 
comparison to V band, we do not use / band data in our 
analysis. However, I band phase light curves for all the 28 
variables along with their time-magnitude photometric data 
in V and / bands, are available in electronic form. In Fig. 7, 
we provide a finding chart of a ~ 13 x 13 arcmin field of the 
cluster marking the variable stars. 



4.3 Comparison with other catalogues and 
previously known variables 

No time-series photometric observations of NGC 6866 have 
been published prior to MOL09. They detected 19 variables, 
12 of which were periodic variables. Out of 19 variable stars, 
we are able to find periodicity in 12 variables while 2 irreg- 
ular variables V15 (ID 53) and V16 (ID 89) in their list has 
also shown irregular variation in our data. Remaining five 
variables from MOL09 are not included in our list of vari- 
able stars. Variable V2 (ID 30) is reported as 5 Scuti star 
with an amplitude variation of only 0.006 mag. We are un- 
able to detect such a low-amplitude variation in this bright 
(V ~ 12.3 mag) star owing to the uncertainties in our data 
towards brighter end. Star V08 (ID 92) is heavily blended 
by two very nearby stars brighter than this star. Though the 
light curve shows a significant flux variation in our data but 
not an ellipsoidal in nature as reported by MOL09. Since 
this is a contaminated star, its variability is highly suscep- 
tible to include it in our list of variable stars. The star V10 



(ID 82) is listed as an eclipsing binary (EA) in their cata- 
logue. For a period of 1.916 d, as reported by MOL09 for 
this star, we should be able to confirm the variations because 
our observations extend over a longer time span. However, 
we cannot confirm the eclipsing nature of this star from the 
present data. Two stars in their catalogue, V17 and V19, 
which were classified as irregular variables lie close to the 
edges of our frames. Since their reliable photometry could 
not be determined in at least 100 frames, we exclude them 
from our list of variable stars. 

It can be seen from the Table 4 that there is a signifi- 
cant difference between our periods and those of MOL09 for 
5 variables. We therefore created phase-folded light curves 
using periods from MOL09 for each variable in common and 
compared them by eye with our phased light curves. We do 
not see any systematic periodic variation for these variables 
using their period apart from the 5 Scuti star V3 which 
could have either period. Variable V4 (ID 1274) is a faint 
star (V ~ 19.1 mag) which lies close to edge of our CCD 
chip in most of the frames hence could be detected only in 
111 frames. We therefore caution that the our estimated pe- 
riod may not be very reliable due to its small number of 
data points. 



4.4 Cluster membership of the variable stars 

We assigned a probability of cluster membership to a star 
on the basis of its angular distance from the cluster center, 
its location in the color-magnitude diagram and its proper 
motion with reference to the mean proper motion of the 
cluster. The spatial probability (P sp ), which gives the posi- 
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Figure 7. Finding chart for the variable stars in the field of NGC 6866. The size of points signifies the magnitude of the stars. 28 periodic 
variables identified in the present study are additionally indicated with open circles. 



tion of the star in the field of the cluster, was determined as 
1 — Ti/r c , where r% is the angular distance of stars from the 
cluster center, and r c (7 arcmin) is the radius of the cluster. 

The photometric probability (P p h) was computed with 
reference to the blue and red limits in the B — V vs V dia- 
gram. The blue sequence was defined using empirical ZAMS 
colors (Schmidt-Kaler 1992) shifted in magnitude and color 
by 11.15 and 0.10 mag respectively to visually match the 
cluster sequence. The red sequence was defined by a shift 
of -0.75 mag in V and a shift of 0.042 in B - V in or- 
der to account for unresolved MS binaries (Maeder 1974, 



Kharchenko 2004). Stars with B — V colors lying within 
the binary sequence of the MS were assigned P p h of 1 and 
are probable members of the cluster, while stars deviat- 
ing along either direction were assigned a probability as 
cxp[-0.5 x (A(B - V)/a {B - V) ) 2 ], where A(B - V) is the 
difference of color from blue or red color limits and <T( S _v) 
is the photometric error in color. 

The proper motion probability (P pm ) was calculated 
using the catalogue of Roeser et al. (2010) which lists po- 
sitions in ICRS system (accuracy 80 to 300-mas) and ab- 
solute proper motion (accuracy 4 to 10-mas/yr) of about 
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Figure 8. Proper motion distribution (left panel) and corresponding CMD (right panel) in the field of NGC 6866. In both panels we 
show stars brighter than 18 mag in V by small dots, while stars used to estimate the mean proper motion are shown by filled circles. A 
circle of 4 mas radius is drawn on the left panel to mark stars having proper motion probability of 1. 



900 million stars, derived from the optical USNOB1.0 and 
near-infrared 2MASS all sky catalogs. A cross-match of our 
photometry (2473 point sources) with that of Roeser et al. 
(2010) gave 1772 stars in common within a radius of 500 
mas. For these stars, a scatter plot of the proper motion in 
right ascension (/x a ) versus that in declination (/«) is plot- 
ted in Fig. 8. Unlike other nearby well studied open clusters 
(see, e.g. Bellini et al. 2010), no clear trend separating mem- 
bers and field stars is seen. We therefore computed the mean 
proper motion of the cluster with V < 17 mag, Pph — 1 
and Psp > 0.5, which gives mean proper motions, U, x = 
-3.47 ± 0.41 mas/yr and p, y = -3.30 ± 0.35 mas/yrQ The 
computations were done iteratively using a er-clipping al- 
gorithm and the uncertainties denote rms deviations. The 
Pph for each star was computed following Kharchenko et 
al (2004) as exp {-0.25 [(jt a - fi a ) 2 /al + (fls - fisf/aj]}, 
where a\ = ct^ q + and o\ = <x?L + o\, . We could thus 
assign P p h for 1949 stars in common with the R10 catalog 
within a search radius of 1 arcsec. The probabilities P ap , P p h 
and P pm are listed in Table 2. 

On the basis of spatial probability, stars found in the 
core of the cluster NGC 6866 (R CO re = 2.0 ± 0.5 arcmin), 
for which P sp is greater than 0.71±0.06, could be the clus- 
ter member. The star having P p h=1.0 might belong to the 
cluster on the basis of photometric probability and, stars 
which have P pm > 60% (l-er) are most probably belong to 
the cluster on the basis of proper motion criteria. Five vari- 
able stars detected in the present study which satisfy all the 

1 The mean proper motion of NGC 6866 is consistent with the 
value (-3.86 ± 0.16, -4.63 ± 0.17) estimated by M09 who calcu- 
lated an older catalogue of Roeser et al (2008) which lists proper 
motion data of only brightest (V < 14 mag) stars. 



three probability criteria are considered as definite members 
of the cluster NGC 6866. Nine stars which satisfy two cri- 
teria are considered as likely members while other 9 stars 
those satisfy either only one criterion or have P pm < 60% 
are considered as unlikely members. Five stars that do not 
follow any membership criteria are considered as definite 
field stars. Membership probabilities and their status as a 
cluster member are given in the Table 4. 

4.5 Characterization of the variable stars 

We assessed the classifications of variable stars by manu- 
ally comparing their phase-folded light curves to template 
light curves of different class of variables. Our classification 
is primarily based on period, shape of the light curves and 
location of the star in the CMD. In Fig. 9, we show the posi- 
tion of the variable stars in the (B — V)o — My plane of the 
cluster. The intrinsic magnitude and color of the variable 
stars are determined using the distance modulus (m — M)v 
= 11.15 mag and extinction E(B — V) = 0.10 mag as es- 
timated in our study. The dashed line is the ZAMS taken 
from Schmidt-Kaler (1982). In order to draw the classical 
instability strip in the CMD, we transformed the theoretical 
instability strip boundaries of Pamyatnysk et al. (2000) into 
our observational plane using color- T e ff relations from Van- 
denBerg & Clem (2003) . The characterization of the variable 
stars detected in our study is summarized as follows: 

4-5.1 Previously identified variables 

ID 27 is classified as a 5 Scuti by MOL09. In the CMD, this 
star is located within the S Scuti instability strip, justify- 
ing their classification. IDs 39 and 94 are reported as 7-Dor 
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Figure 9. The (B - V)o vs M v CMD for NGC 6866. The thick dashed line shown by blue line is the ZAMS corresponding to V - M v 
= 11.15 and E(B-V) = 0.10 mag, while the red dashed line indicates a shift in magnitude and colour due to binarity. Positions of 28 
periodic variables are also shown in the CMD. The standard errors in colour and magnitude at various magnitude ranges are given in 
the right corner. The thin solid lines represent the instability strip for <5-Sct stars while the thick solid lines represent the same for 7-dor 
stars. 



variables. The location of ID 39 is towards the blue edge of 
the S Scuti instability strip, but its period suggests that it 
is a 7-Dor star. This makes it an interesting target for the 
future observations. Star 81 is classified as irregular variable 
in MOL09. We found a periodicity of 1.24 day in the present 
data, however, our observations do not give a uniform sam- 
pling at all the phases hence period may be little uncertain. 
This star along with the star 94 fall in the 7-Dor instability 
strip in the HR diagram which suggests that these stars be- 
long to 7-Dor class of variables. MOL09 classified stars 74, 
487, 494 and 1274 as W Ursa Majoris (W UMa) type vari- 
ables. Our light curves support their classification for the 
stars 487, 494 and 1274, however, light curve of the star 74 
looks like an ellipsoidal variation. Star 248 has an unusual 
light curve. This star is reported as an ellipsoidal variable 
but it does not have the typical light curve of interacting 
binaries. It shows the light minima of equal depth instead 
of the maxima. Such a strange light curve might arise if we 
view an interacting binary in a highly-eccentric orbit close to 
pole-on. The variation in brightness is a result of the change 
in sizes of the stars as they start approach periastron. Stars 
239 and 253 are classified as irregular variables by MOL09, 



but our phased light curves show clear periodic variability. 
Star 239, which is an unlikely member of the cluster, is red- 
dest (B -V = 2.09, V — I = 3.44) and longest period (~ 
37 day) variable identified in the present study. This star is 
most-likely a background red giant and variability in these 
stars shows considerable periodicity in their light variation, 
accompanied or sometimes interrupted by various irregulari- 
ties as seen in its light curve. We therefore classified this star 
as a semi-regular variable. Star 253 seems to be a rotating 
variable. 

Star 47 is reported as 0.066 day <5-Sct variable (VI) by 
MOL09. Though this star falls in the <5-Sct instability strip 
in the CMD of our data, the shape of light curve is not a 
typical 5-Sct variation. Moreover, our study finds a period of 
0.275 day and magnitude variation of 0.041 mag for this star 
which is significantly higher in comparison of MOL09. While 
proper motion probability for this star is reported as 72% 
by MOL09, present study gives a probability of only 55%. 
Based on the spatial distribution, this star has only 27% 
probability of belonging to the cluster. Therefore, despite 
this star falling in the 5-Sct instability strip in the CMD 
of the cluster, present observations suggests that the star is 
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unlikely a cluster member and highly susceptible for being 
a 5-Sct variable. 



4-5.2 New variables 

A total of 16 new variables were detected from the time- 
series data analysis of NGC 6866. We classify star 36, a likely 
member of the cluster, as S Scuti variable as it is located in 
the 5 Scuti instability strip. Star 58, which is a field star, 
shows a light variation similar to a 7-Dor variable. Stars 16, 
20 and 191 show similar light curves. They could be binary 
stars, though further data is needed to ascertain their ex- 
act nature of variability. Star 158 seems to be an eclipsing 
binary of the Algol (EA) type, where the primary and sec- 
ondary eclipses are clearly seen in the phased light curve. 
Stars 221 and 231 show similar light curves. The periods 
and amplitudes are also similar. Star 221 is a likely mem- 
ber, 231 is a confirmed member of the cluster based on their 
membership probabilities. Both stars lie on or near the main 
sequence. These stars are placed in the group of pulsating 
variables (PV). Stars 332 and 349 look like a semi-regular 
or long period variable (LPV) . Star 444 has a periodicity of 
~ 16.3 days and display signature of primary and secondary 
eclipses. More observations will be needed to see light varia- 
tions at all the phases for this long period variable, neverthe- 
less, we classify this star as /3 Lyr (EB) type binary based 
on present observations. Stars 1077, 1088, 1292, 1421 and 
1583 have periods similar to S Scuti stars, but have signif- 
icantly large amplitudes, suggesting that these stars could 
be high-amplitude <5-Sct stars (HADS). 

Present study along with MOL09 classified stars 487, 
494 and 1274 as W UMa type variable which are known 
to obey a PLC relationship between the period, color and 
absolute brightness (Rucinski 2004) as follows: 

M v = -4.44Zop(P)+3.02(B-U) +0.12 a = 0.25 mag 

where (B — V)o is the intrinsic color index and P is the or- 
bital period in days. We derived the apparent distance mod- 
ulus (m — M)v of each W UMa star from the absolute mag- 
nitude Mv determined using the above relation and their 
mean V magnitudes. A reddening of E(B — V) = 0.10 was 
assumed. Of the three W UMa variables, star 487 is the most 
interesting. It has a very high proper motion probability of 
0.85, spatial probability of 0.76 and photometric probability 
of 1.0 and hence considered as a confirmed member of the 
cluster. However, (m — M)v of this star is estimated as 12.53 
mag from the W UMa PLC relationship. This is 1.38 mag 
fainter than the apparent distance modulus of the cluster. 
Since the period and colour of this star is quite robust, a 
large value of (m — M)v suggests that either this star is not 
a cluster member or it does not belong to class of W UMa 
variables hence not following the above relationship. This 
star needs further attention to ascertain its true member- 
ship and nature of variability. Other two stars 494 and 1274 
has zero proper motion and no proper motion information 
available respectively. Though both of them fall in the CMD, 
their spatial position in the target field is quite far from the 
cluster center. Furthermore, stars 494 and 1274 are fainter 
by 0.94 and 2.03 mag respectively than the estimated dis- 
tance modulus of the cluster. This makes these two W UMa 
stars unlikely members of the cluster. 



5 SUMMARY 

The work presented here is the first paper of our series of 
papers on detection and characterization of variable stars 
in young and intermediate age open clusters. In this work 
we have presented a search for variable stars in the cluster 
NGC 6866. We found 28 variables in the period range 0.03- 
day (~ 48-mins) to 37-day and confirm the irregularity in 
the light curves of other two irregular variables reported 
by MOL09. Among them, 16 are newly-discovered periodic 
variables. We also determined period for 3 other variables 
which were reported as irregulars in MOL09. Prior to this 
study, the shortest period variable detected in NGC 6866 
was a 13 mag <5-Sct star with a period of 1.6-hr. Since we 
have carried out continuous observations in the V band for 
more than 4 hours on 3 nights, this has enabled us to find 
variable star with period as short as 48 min. In the field of 
NGC 6866, we did not find any variable star which shows 
variations in excess of 0.5 mag. In comparison with previous 
studies, we report more than twice the number of periodic 
variables in the field of NGC 6866. 

In the present study, we have analysed membership of 
the stars on the basis of their distance from the cluster cen- 
ter, positions on the CMD, and their proper motions, wher- 
ever available. The membership probabilities of the stars 
on the basis of all the three criteria were estimated. Our 
analysis of 28 periodic variables suggests that 14 variables 
are either confirmed members or likely members of the clus- 
ter, while remaining 14 variables are either field stars or 
unlikely members of the cluster. Based on the shape and 
period of the light variations, together with the colour and 
amplitude, we found several 5 Scuti, 7-Dor, rotational vari- 
able and eclipsing binaries. Few variables could not be clas- 
sified with precision and more photometric observations are 
needed to ascertain the exact nature of these stars. Further 
multi-colour photometry and spectroscopic observations of 
these stars will help to determine their parameters such as 
masses and radii. 

We also provide calibrated UBVRI photometry of 2473 
stars down to V ~ 21.5 mag. These data were used to 
determine the physical parameters of the cluster and de- 
rive the cluster radius from the stellar density profile and 
the extinction, distance and age from the colour-colour and 
colour-magnitude diagrams. The basic parameters of the 
cluster NGC 6866 is obtained through isochrone fitting, giv- 
ing log(t) = 8.8 yr, a distance modulus of (m — M)o = 10.84 
mag, and extinction E(B — V) = 0.10 mag. The radial distri- 
bution of the stellar surface density indicates that the core 
and cluster radius is extended up to about 2 and 7 arcmin 
respectively with a peak density of 5.7±0.7 star/arcmin 2 at 
the cluster center. 
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